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Uniformly rotating axisymmetric fluid conflgurations 
bifurcating from highly flattened Maclaurin spheroids 
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ABSTRACT 

We give a thorough investigation of sequences of uniformly rotating, homogeneous 
axisymmetric Newtonian equihbrium configurations that bifurcate from highly flat- 
tened Maclaurin spheroids. Each one of these sequences possesses a mass-shedding 
limit. Starting at this point, the sequences proceed towards the Maclaurin sequence 
and beyond. The first sequence leads to the well known Dyson rings, whereas the end 
points of the higher sequences are characterized by the formation of a two-body sys- 
tem, either a core-ring system (for the second, the fourth etc. sequence) or a two-ring 
system (for the third, the fifth etc. sequence). Although the general qualitative picture 
drawn by Eriguchi and Hachisu in the eighties has been confirmed, slight differences 
turned out in the interpretation of the origin of the first two-ring sequence and in the 
general appearance of fluid bodies belonging to higher sequences. 
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1 INTRODUCTION 



If one moves along the Maclaurin sequence of uniformly 
rotating, ajcisymmetric, homogeneous fluid ellipsoids with 
fixed mass and fixed density, starting at the non-rotating 
configuration and proceeding towards increasing angu- 
lar momentum, one first encounters the Jacobi branch 
point where the ellipsoids become secularly unstable with 
regard to the first non-axisymmetric perturbation (see 
e.g. Chandrasekhar 1969). The corresponding Jacobi se- 
quence branching off at this point leads to further bifur- 
cations, in particular to Poincare's 'pear-shaped' sequence 
(see Eriguchi & Hachisu 1982b). Farther along the Maclau- 
rin path, one next comes to the bifurcation points of the 
non-axisymmetric 'triangle', 'square' and 'ammonite' se- 
quences (Chandrasekhar 1969, Eriguchi & Hachisu 1982a), 
before one arrives at the first axisymmetric sequence bifur- 
cating at an eccentricity ei = 0.98523 (Chandrasekhar 1967, 
Bardeen 1971). As conjectured by Bardeen (1971) and con- 
firmed by Eriguchi & Sugimoto (1981), the bodies of this 
sequence pinch together gradually at the centre and eventu- 
ally form the anchor-ring configurations studied by Dyson 
(1892, 1893) and Wong (1974); see also Poincare (1885), 
Kowalewsky (1895), Lichtenstein (1933) and a somewhat re- 
lated paper by Kley (1996). 

Another not yet confirmed conjecture by Bardeen 



E-mail:ansorg@tpi.uni-jena.de 



(1971) concerns a sequence of axisymmetric 'central-bulge' 
configurations hkewise bifurcating at £1. In this paper we in- 
deed found this sequence, which might be considered to be a 
continuation of the Dyson-ring sequence beyond the Maclau- 
rin ellipsoids and ends in a mass-shedding limit. However, 
the surface shapes of the corresponding fluid bodies do not 
generate a 'central-bulge' region; their appearance is more 
'lens shaped'. 

Apart from the study of the Dyson-ring sequence we 
shall give a detailed analysis of the next axisymmetric se- 
quences bifurcating from the Maclaurin ellipsoids. The first 
core-ring sequence branches off at 62 ~ 0.99375. Starting at 
a mass-shedding limit, this sequence proceeds towards the 
Maclaurin sequence and beyond, finally leading to the for- 
mation of a core-ring system. 

In contrast to the results by Eriguchi & Hachisu (1982a) 
who stated that the first two-ring sequence branches off at 
£2, we found that the bifurcation occurs at £3 = 0.99657. 
Again, a mass-shedding limit marks one end point of this 
sequence, leading from here towards and beyond the Maclau- 
rin sequence and ending in the formation of a two-ring sys- 
tem. 

The same qualitative picture repeats as one moves 
to the higher axisymmetric bifurcation points at £4 = 
0.99784, £5 = 0.99851 etc., of which there are infinitely 
many, accumulating at £ = 1 (Bardeen 1971). Starting at 
a mass-shedding limit, the sequences proceed towards the 
Maclaurin sequence and beyond, leading eventually to the 
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Table 1: Physical quantities for the Maclaurin ellipsoids at the 



bifurcation points Sk- 
Section 4. 



For the definition of i 



T and W see 



k 




Axis ratio 


a;2 


f 


T/\W\ 


1 


0.98523 


0.17126 


0.087262 


0.021741 


0.35890 


2 


0.99375 


0.11160 


0.066105 


0.029152 


0.40345 


3 


0.99657 


0.082750 


0.052711 


0.034638 


0.42664 


4 


0.99784 


0.065744 


0.043714 


0.039037 


0.44084 


5 


0.99851 


0.054534 


0.037301 


0.042740 


0.45044 


6 


0.99891 


0.046589 


0.032513 


0.045957 


0.45736 


7 


0.99917 


0.040664 


0.028806 


0.048814 


0.46259 


8 


0.99935 


0.036075 


0.025854 


0.051395 


0.46667 


9 


0.99947 


0.032417 


0.023449 


0.053755 


0.46995 


10 


0.99957 


0.029432 


0.021451 


0.055935 


0.47265 



formation of a core-ring system for 621 and a two-ring system 
for 621+1, see Figs 7-11. The body's surface is characterized 
by a particular number of grooves, notably I, for configura- 
tions close to the two-body systems that are the end-stages 
of the sequences bifurcating at 821 and 621+1- The outer- 
most groove pinches together first and an outer ring without 
grooves separates. In some continuation process beyond this 
formation of a two-body system, the other grooves might 
also pinch off, leading eventually to a multi-body system, at 
most consisting of a core and I rings or {I + 1) rings respec- 
tively]^ However, since a continuation process of this kind is 
not unique, we consider the sequences in question precisely 
up to the formation of the two-body system. 

These investigations have been carried out up to the 
eio-sequence. The corresponding final configurations, when 
the two-body system is formed, can be seen in Figs 7-11. 

Representative plots of physical quantities as well as of 
meridional cross-sections have been provided up to the £5- 
sequence, see Figs 1-10. Note that the Maclaurin sequence 
becomes dynamically unstable with respect to axisymmetric 
perturbations for e > 0.99856 (Bardeen 1971, Eriguchi & 
Hachisu 1985). Therefore, if we restrict our considerations 
to axial symmetry, the sequences bifurcating at ei, £2, ■ • ■ , £5 
are more relevant than those branching off at Eg, £7, . . . {sk ^ 
£6 = 0.99891 for > 5). 

Table 1 lists physical quantities for the Maclaurin ellip- 
soids at the bifurcation points Ek for fc = 1, 2, . . . , 10. Tables 
2-4 contain numerical data with an accuracy of five digits 
for configurations of the £1, £2 and £3-sequences. Addition- 
ally, Table 5 lists the ratios of the masses of the inner to the 
outer body at the two-body formation point for the £2 . . . £5- 
sequences. 



2 BIFURCATION POINTS OF 
AXISYMMETRIC SEQUENCES 

In this Section we review the results by Bardeen (1971) 
who obtained all bifurcation points on the Maclaurin se- 
quence which correspond to axisymmetric perturbations. 
For our purposes, we will use the Maclaurin ellipsoids at 



^ This result is in contrast to the statement by Eriguchi & 
Hachisu (1982a), who expected the formation of k rings for the 
Efc-sequence. 



these points as initial solutions for finding the respective 
sequence branching off. 

The solution <1> of the axisymmetric Poisson equation 



(1) 



where G is Newton's constant of gravity, fi the mass density, 
and A the Laplace operator, given in cylindrical coordinates 
(P,C) by 

A- — + -— + — 

dp^ p dp dC,^ 

can be written in the form (Morse & Feshbach 1953) 

00 

<E-(e, 77) = -2™g G ^(2; + 1)P, (r,) X 

i=0 



i -1 

e 1 



Here, the oblate spheroidal coordinates ^ and 77 are con- 
nected with p and C, by 



p = aoV(l + e)(l-'72) , C = ao^^- 

The functions Pi are Legendre polynomials, and the func- 
tions pi and qi are defined by 



(i+c'^)h(^')P 



The fluid body revolving with the uniform angular velocity 
Q is in hydrostatic equilibrium if the condition 



^iOv) - \^^P^ = $0 = const. 



(2) 



is satisfied at every point of the body's surface. 

In order to derive the Maclaurin solutions, one simply 
starts with an ellipsoidal shape of the surface characterized 
by (,s{ri) = Co = const, and moreover takes ao to be the 
corresponding focal length. The density p vanishes for f > 
and is a positive constant if ^ < .^o which in the following 
we simply call p. Then the resulting surface potential reads 
as follows 

*(?o,»7) = -^™oGMeo(l + Co)[go(?o)+g2(Co)P2(r;)], 

and the surface condition (^) indeed is satisfied, provided 
that Q. and $0 are well-defined function of ao, and p which 
can be read off from (|^). 

If we now consider an axisymmetric perturbation of the 
above Maclaurin shape 

Cs(7)) = eo + <i(77) + o('^'), 

insert the corresponding surface potential into the condi- 
tion and linearize the resulting expression with respect 
to the small parameter k, we uniquely find a neighbouring 
Maclaurin ellipsoid unless 



2p2fc+2(Co)l72fc-|-2(Co) 



Co 



P2(Co) 



[1-2^092 (Co)] 



(3) 
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Figure 1. For the first five axisymmetric sequences (1), . . . , (5) , the squared angular velocity up" = n^/(47rG/i) is plotted against the 
radius ratio ri/p2, where ri = (^i ((^i : polar radius) for spheroidal figures and ri = — pi (pi : inner equatorial radius) for toroidal shapes; 
P2 is always the radius of the outer equatorial rim. The dotted curve represents the Maclaurin sequence and the dashed one corresponds 
to the Dyson approximation. 
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Figure 2. For the first five axisymmetric sequences, the squared angular velocity = / (AnGfj.) is plotted against the dimensionless 
squared angular momentum j^, given in formula (pi). Dotted and dashed curves again refer to the Maclaurin sequence and the Dyson 
approximation respectively. The •'s mark the bifurcation points on the Maclaurin sequence and □ the transition configuration of spheroidal 
to toroidal bodies on the Dyson-ring sequence. 
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is valid for some fc£{l,2,3,...}. The corresponding eccen- 
tricity £fc belonging to the root ^q''' of the above equation is 
given by 



k = 1,2,3, 



and marks the bifurcation point on the Maclaurin sequence 
at which the fc**^ axisymmetric sequence branches off. The 
numerical values of the first ten Sk are given in Table 1. 
As illustrated in Section 4, the bifurcation points appear as 
cross-over points of the respective sequence with the Maclau- 
rin sequence if one plots the evolution of these sequences in 
an appropriate parameter diagram (such as angular velocity 
against the ratio of polar and equatorial radii). Note that 

(k) 

the corresponding linear surface displacement ^{ turns out 
to be 

P2k+2{ll) 



(4) 



3 NUMERICAL METHODS 

The Poisson-integral 



$(r) = -Gfi 
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over the fiuid body V can be reduced to a one-dimensional 
integral as foUows (HiU & Wheeler 1953, Wong 1974): 

^(P.C) 

P2 Cs(p') 27r 

= -Gil 

PI -Cs(p') 

P2 Cs{p') 

= -4Gfi 



p'd<f)'dC'dp 



yjp^ + P'^ + (C - C? - 2pp'cos0' 
p'K{m)dCdp' 



PI -Cs(p') 

P2 



(5) 
(6) 



= -2Gp J [A-{p';p,0 - A+{p';pX)]p'dp' 
PI 

with 

A±{p';p,0 = 

±[{p + p')K{m±) ~ 2pD{m±)]C{p') + [C T Up'^jm^) 

In these expressions the surface is described by a non- 
negative function ("a defined on the interval [pi , p2] and ("^ 
denotes its first derivative. Note that pi — for a spheroidal 
equlibrium figure, whereas for a toroidal figure pi denotes 
the inner equatorial radius; p2 is always the radius of the 
outer equatorial rim. 

The functions K and D are defined by the elliptic inte- 
grals 



7r/2 



K{m) 



da 



m sin a 



Figure 3. Magnified details of Fig. 1. The axes of the abscissae 
and ordinates have been stretched by a diff'erent factor. 



Tr/2 



D{m) 



sin a da 



\/ 1 — sin^ a 



and the modulus is given by 



4pp' 



and 



m± 



App' 



[p + P? + [(TCsip'W 
in the expressions ^ and respectively^ 



^ Writing as a contour integral over a vector field and using 
Stoke's integral theorem as well as functional relations between 
elliptic integrals, one arrives at the expression ra). 
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Figure 4. Magnified details of Fig. 2. The axes of the abscissae 
and ordinates have been stretched by a diff'erent factor. 



Inserting the integral expression (^) into the surface 
condition (0), one gets a non-linear integral equation for 
the unknown shape of the body's surface, which is uniquely 
soluble at least in the vicinity of some given solution if an 
appropriate triple of parameters is prescribed^. We write 
the shape of the body's surface parametrically in terms of a 
Chebyshev expansion in the following manner: 



Ps{r) 

9{r) 
Tk{x) 



y/p'i + (P2 - Pl)r , 



V(l-r)[C?+rff(r)], 

k = l 

cos(fc arccos x). 



The variable r runs over the interval [0, 1], and the parame- 
ter (i represents the polar radius for spheroidal figures and 
(^1=0 for toroidal shapes. 

If we prescribe /i, Q, and $o, restrict ourselves to a finite 
number (A'' — 2) of coefficients Qk and take the surface con- 
dition (Q) at the A'' surface points (ps(Tfc), Csfps (rfc)]) with 
TN{2Tk — 1) = 0, the integral equation in question turns 
into a system of A'^ equations for the unknown coefficients 
gk and the likewise unknown quantities ri and p2 where ri 
stands for (^i in the spheroidal and for (—pi) in the toroidal 
case. Alternatively, we can prescribe the value g{T — 1) and 

As an example, the density fi, the angular velocity fl and the 
constant of the fluid body can be taken. However, at the bi- 
furcation points Ek it is necessary to use different parameters, in 
order to assure a unique solution, see below. 
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Figure 5. For the Dyson-ring sequence, the dimensionless total 
energy E, given in formula is plotted against 47rj^, see (|^). 
The Maclaurin sequence is indicated by the dotted curve; the 
• marks the corresponding bifurcation point. The diagram (b) 
shows a magnified detail of (a), again with different stretching 
factors for the two axes. In (c), the difference (E — Ej^ji) of the 
respective energies of Dyson-ring configurations and Maclaurin 
solutions belonging to the same is displayed. The □ marks the 
transition configuration of spheroidal to toroidal bodies. 



search instead for $o. A choice g{l) 7^ assures a unique 
solution at the desired £fe-sequence if we are in the vicinity 
of Sk- This is due to the fact, that for all Maclaurin solutions 
gij) is identically zero. 

The solution of the remaining system of TV equations 
can be found by means of a Newton-Raphson method. The 
Maclaurin ellipsoids at the bifurcation points serve as ini- 
tial solutions, and by choosing g{l) < we move towards 
the mass-shedding limit of the respective sequence whereas 
g{l) > yields the path towards the ring or core-ring struc- 
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Table 2: Physical quantities for the ei -configurations (A) . . . (L) 
displayed in Fig. 6. The abbreviations 'MSL' and 'MLE' refer to 
the mass-shedding limit and the Maclaurin ellipsoid respectively. 



t\Ipi 






TI\W\ 






0.19224 


0.092400 


0.019847 


0.34379 


(A): 


MSL 


0.18500 


0.089834 


0.020815 


0.35175 


(S) 




0.17126 


0.087262 


0.021741 


0.35890 


(C): 


MLE 


0.16000 


0.085614 


0.022260 


0.36249 


(O) 




0.13000 


0.082464 


0.022927 


0.36549 


{E) 




0.10000 


0.080648 


0.022898 


0.36260 


(F) 




0.050000 


0.079809 


0.022215 


0.35404 


(G) 




0.000000 


0.080181 


0.021836 


0.35047 


(H) 




-0.20000 


0.078338 


0.020036 


0.32825 


(I) 




-0.40000 


0.057176 


0.021726 


0.31082 


(J) 




-0.60000 


0.030207 


0.031052 


0.31542 


(K) 




-0.80000 


0.0089862 


0.063778 


0.33800 


(L) 





tures, see Section 4. The mass-shedding limit is reached 
when g{l) = — Ci, which leads to a cusp at the equatorial 
rim of the body. Although for this configuration derivatives 
of the above function g are singular at r = 1 (which results 
in a loss of accuracy, see below), we can choose (/(1)/Ci a-s 
one of our parameters and can thus place ourselves precisely 
on this point. 

Once upon the sequence and sufficiently far away from 
the Maclaurin sequence, we may change the parameter pre- 
scription, e.g. to the triples (/i, Q,, "I>o) or (/i, ri/p2, ^o)- Note 
that rapid convergence of this method towards the actual 
shape of the body as A'' grows, requires, for a given set of 
coefficients ot, a highly accurate evaluation of the integral 
expression (Bl). This is done by dividing the interval [0, 1] into 
many subintervals, which are extremely small in the vicin- 
ity of the respective zero Tk and grow exponentially as one 
moves away from Tk- Moreover, for m ~ 1 the elliptic func- 
tions are written in logarithmic expansions (see Gradstein & 
Ryshik 1981, p. 295 f.). Additionally, special care is taken for 
surface points close to the equatorial plane. In this manner 
we achieve an accuracy of up to 12 digits for intermediate so- 
lutions that are sufficiently far away from the mass-shedding 
limit as well as from any transition body where the topology 
of the figure changes. The loss of accuracy as one approaches 
these critical configurations comes about since derivatives of 
g become singular at a certain point^on the body's surface. 
However, we are still able to obtain an accuracy of 5 digits 
for all physical quantities given in Tables 2-5. 



This critical point is the equatorial rim for the mass-shedding 
configuration, the 'pinching' point for the configuration at which 
the two-body system forms and the origin p = 
transition body of spheroidal to toroidal figures. 



Table 3: Physical quantities for the e2-configurations (A) . . . (L) 
displayed in Fig. 7. The abbreviations 'MSL', 'MLE' and '2-BS' 
refer to the mass-shedding limit, the Maclaurin ellipsoid and the 
formation of the two-body system respectively. 



ri/p2 






T/\W\ 






0.11207 


0.069411 


0.027909 


0.39710 


{A): 


MSL 


0.11140 


0.068335 


0.028310 


0.39919 


(B) 




0.11160 


0.066105 


0.029152 


0.40345 


{cy. 


MLE 


0.11300 


0.064738 


0.029663 


0.40587 






0.11500 


0.063529 


0.030101 


0.40779 


(E) 




0.12000 


0.061475 


0.030787 


0.41032 


(n 




0.12500 


0.059960 


0.031218 


0.41132 


(G) 




0.13000 


0.058718 


0.031496 


0.41137 


{H) 




0.13500 


0.057661 


0.031659 


0.41067 


(I) 




0.14500 


0.055961 


0.031701 


0.40744 


(J) 




0.15500 


0.054729 


0.031417 


0.40219 


(K) 




0.17558 


0.054327 


0.029720 


0.38610 


(L): 


2-BS 



Table 4: Physical quantities for the es-configurations (A) . . . (L) 
displayed in Fig. 8. The abbreviations 'MSL', 'MLE' and '2-BS' 
refer to the mass-shedding limit, the Maclaurin ellipsoid and the 
formation of the two-body system respectively. 



ri/p2 




f 


T/\W\ 






0.087646 


0.054735 


0.033738 


0.42327 


(A): 


MSL 


0.086000 


0.053836 


0.034134 


0.42477 


(B) 




0.082750 


0.052711 


0.034638 


0.42664 


(C): 


MLE 


0.075000 


0.050871 


0.035462 


0.42947 


(D) 




0.060000 


0.048799 


0.036307 


0.43162 


(E) 




0.040000 


0.047554 


0.036625 


0.43116 


(F) 




0.020000 


0.047289 


0.036544 


0.42979 


(G) 




0.000000 


0.047391 


0.036458 


0.42932 


(H) 




-0.050000 


0.047666 


0.036216 


0.42800 


(I) 




-0.10000 


0.046297 


0.036022 


0.42281 


(J) 




-0.12000 


0.045352 


0.035860 


0.41895 


(K) 




-0.13998 


0.044952 


0.035174 


0.41220 


(L): 


2-BS 



4 RESULTS 

The evolutions of the five sequences bifurcating at the ec- 
centricities El, . . . , £5 are displayed in the Figs 1 and 2 where 
oj^ — f2^/(47rG/i) is plotted against the ratio r\/p2 and 
against the dimensionless angular momentum defined by 

f = -La (7) 

where M and J are the mass and the angular momentum 
of the body, respectively. The corresponding (dotted) curves 
of the Maclaurin ellipsoids are given as well as the analytic 
Dyson approximation of the ei-sequence (dashed lines). 

In Fig. 1, the bifurcation points appear as cross-over 
points of the Maclaurin and respective Sfc-sequence. One end 
of a particular sequence is marked by a mass-shedding limit, 
always possessing a higher angular velocity lo^ compared 
to the other end. The details of the typical evolution of a 
particular sequence can be seen in Fig. 3. A characteristic 
feature of all sequences is the occurence of a local maximum 
in close to the mass-shedding limit, see Figs 3(b), 3(c) 
and 4(b). 

= C for the AH £2!-i-sequences show the transition of a spheroidal 

to a toroidal figure, and u)'^ possesses a local minimum and 
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Figure 6. Meridional cross-sections of configurations belonging 
to the ei-sequence, see Table 2. The dimensionless quantity C/P2 
is plotted against p/p2- The Maclaurin body is denoted by hatch- 
ing. 




Figure 7. Meridional cross-sections of configurations belonging 
to the £2-sequence, see Table 3. The dimensionless quantity C/P2 
is plotted against p/p2- The Maclaurin body is denoted by hatch- 
ing. 



maximum before and after this point respectively, see also 
Figs 3(a) and 4(a). Another typical property is the vanishing 
gradient of lo'^ as a function of ri / p2 as well as a discontin- 
uous change of the sign of its curvature at r\/ p2 = 0. Apart 
from the £i-sequence, the paths end with the formation of 
a two-body system, namely a core and a ring for k = 21 
and two rings if fc = 2/ -I- 1. Shortly before coming to this 
end, a minimum in lj^ occurs, which becomes less and less 
pronounced as one moves to higher values of k. 



As depicted by Hachisu & Eriguchi (1983), the paths 
of the Efe-sequences glance off the Maclaurin sequence^, if 
one plots Lu^ against , see Fig. 2. Representative details 
corresponding to the ei-curve can be seen in Fig. 4. Note the 
extraordinary accuracy needed in order to resolve the global 
maximum in lu'^ close to the mass-shedding limit displayed 
in Fig. 4(b). 



^ At the bifurcation points the slopes are identical. 
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(C) 



(E) 



(F) 



(G) 



(H) 



(I) 



(J) 



(K) 




Figure 8. Meridional cross-sections of configurations belonging 
to the £3-sequence, see Table 4. The dimensionless quantity C/P2 
is plotted against p/p2 - The Maclaurin body is denoted by hatch- 
ing. 



Figure 9. Meridional cross-sections of configurations belonging 
to the e4-sequence. The dimensionless quantity C/P2 is plotted 
against p/p2- The Maclaurin body is denoted by hatching. 



In Figs 1 and 2, moreover the analytic curves of the 
Dyson approximation can be found. With the assumption 
of circular cross-sections of the toroids as p\/p2 1, one 
mav expand ^ in terms of (1 — pij p2) and find by means 
of (y) an expansion of the toroids' surfaces. Dyson went to 
the order o[(l — P1/P2)''], resulting in excellent agreement 
with the numerical result^ for large pi/p2- 

Fig. 5(a) is a replot of fig. 1 from Bardeen (1971), with 
the Dyson approximation being replaced by the actual evolu- 
tion of the £i-sequence. Here, the dimensionless total energy 
E is given by 



E = TT^ TTT = 2-K'IUJ (1 H ) 

with the kinetic and gravitational energies 

r = \.m and W = ^ I <^ d^v = ^M$o + T. 



(8) 



which is valid for all stationary rotating fluid bodies, and 
from the hydrodynamical Euler equation giving for a uni- 
formly rotating homogeneous configuration the pressure p 
by 



m(*o 



Fig. 5(a) led Bardeen to the conjecture that the Dyson- 
ring sequence branches off from the Maclaurin sequence ex- 
actly at the El-bifurcation point, marked in Fig. 5 by •. 
However, a magnified detail of the plot in question shows 
a completely different picture, as can be seen in Figs 5(b) 
and 5(c). Coming from the toroidal end, the ei-sequence 
intersects the Maclaurin sequence close to but not exactly 
at the bifurcation point. As can be seen in Fig. 5(c), the 
ei-sequence intersects the Maclaurin sequence another time 
before it moves on to the bifurcation point and beyond to the 
mass-shedding limil|^ It is rather surprising that from this 



The latter formula follows from the virial theorem 



2T + W + 3 / pd^r = 



^ In our computations, we compared configurations with the 
same pi,P2 and p. 



^ The relation dE = inuj dj , which is valid along any equilibrium 
sequence of uniformly rotating homogeneous fluid bodies (cf . Har- 
tle &; Sharp 1967, Bardeen 1971), explains two features of Fig. 5: 
the two cusps and the fact that the slope of the si-sequence is 
identical to that of the Maclaurin sequence at the bifurcation 
point. Moreover, together with the coinciding slopes in the oj^- 
j'^-plane (Fig. 2), one can conclude that the second derivatives in 
the B-j'^-plane (Fig. 5) are identical at the bifurcation point as 
well. 
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Table 5 : The ratios M\ /Af 2 of the masses of the inner to the outer 
body at the two-body formation point for the £2 ■ • • es-sequences. 



Figure 10. Meridional cross-sections of configurations belonging 
to the £5-sequence. The dimensionless quantity C/P2 is plotted 
against p/p2- The Maclaurin body is denoted by hatching. 




(2) 



(3) 



(4) 



(5) 



M1/M2 0.41051 0.86563 1.7235 2.4987 



Kleinwachter & Meinel (2002). Additionally, the values of 
physical quantities like and are similar for the two 
bodies^. Such a neighbourhood can also be found for the 
configurations [F) and fig. 1 (b) of Ansorg et al. (2002). 

Likewise, Tables 3 and 4 list data for the fiuid bodies 
shown in Figs 7 and 8. Again, there is a configuration (B) 
of the e2-sequence (see Fig. 7) which is in the vicinity of a 
relativistic body, (c) in fig. 1 of Ansorg et al. (2002). Note 
the hint of a minimum at the north pole of the relativistic 
figure. Also for the configuration (B), this property is sug- 
gested by the linear surface displacement (U|). But, as can 
be seen in Figs 6-10, the Newtonian configurations in be- 
tween the Maclaurin body and the mass-shedding limit nei- 
ther generate a central bulge (which one would expect for 
the e2!_i-sequences) nor a central minimum region (which 
would follow for the e2rsequences). However, in the gen- 
eral relativistic vicinity of these solutions, at least the latter 
property seems to show up. 

Table 5 lists the ratios of the masses of the inner to the 
outer body at the two-body formation point for the 62 ■ ■ ■ £5- 
sequences. 

Figs 9 and 10 display configurations of the £4 and £5- 
sequences, and Fig. 11 is devoted to the formation of the 
final two-body-system for the . . . £io-sequences. As men- 
tioned in the Introduction, for £2; the resulting configuration 
consists of a corrugated central core-region and an outer ring 
which just pinches off. Likewise, for £2i+i a two-ring system 
is formed. 

An interesting open question concerns the behaviour of 
the £fc-sequences as k tends to infinity, corresponding to the 
disk limit (e 1) of the Maclaurin ellipsoids. 
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Figure 11. The formation of the two-body system at the end of 
the £6 • • ■ £io-sequences. 



entangled picture the correct connections of the Dyson-ring 
and the Maclaurin sequence could be drawn. 

Table 2 contains numerical data, accurate to five dig- 
its, for physical quantities belonging to the £i-sequence. The 
corresponding meridional cross-sections of these bodies are 
displayed in Fig. 6. The intermediate Maclaurin body (C) 
is denoted by hatching. Note that a lens-shaped configura- 
tion (B) appears that looks similar to the relativistic ho- 
mogeneous fluid body displayed in flg. 1 (a) of Ansorg, 
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